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Abstract 

The suitability of various implementations of inertial-electrostatic confinement (IEC) sys- 
tems for use as D-T, D-D, D- 3 He, p- u B, and p- 6 Li reactors is examined; these IEC designs 
create a deep electrostatic potential well within the plasma in order to confine and acceler- 
ate ions, and they typically use magnetic fields or electrostatic grids to confine electrons. It 
is shown that while an IEC reactor would have the advantages of high power densities and 
relatively simple engineering design when compared with other fusion schemes, it suffers 
from several flaws. Foremost among these problems is that on the basis of velocity-space 
diffusion calculations, it does not appear to be possible for the dense central region of a 
reactor-grade device to maintain significantly non- Maxwellian ion distributions or to keep 
two different ion species at significantly different temperatures; this discovery contradicts 
earlier claims of the particular suitability of IEC systems for advanced fuels. Since the 
ions form a Maxwellian distribution with a mean energy not very much smaller than the 
electrostatic well depth, ions in the energetic tail of the distribution will likely be lost at 
rates greatly in excess of the fusion rate. Furthermore, even optimistic assumptions about 
the performance of the surrounding polyhedral cusp magnetic field lead to the conclusion 
that the electron losses from the machine will be intolerable for all fuels except DT. If 
electrostatic grids are used instead of magnetic cusps, the particle losses should be even 
worse. Finally, because the Maxwellian ion distributions are not fundamentally different 
than those in other fusion schemes, bremsstrahlung losses will be similar; in particular, 
the ratio of bremsstrahlung power to fusion power will at best be ~1.7 for p- 11 B and 
~5.4 for p- 6 Li, demonstrating that IEC cannot utilize proton-based fuels in a practical 
manner. In order for IEC systems to be used as fusion reactors, it will be necessary to 
find methods to circumvent these problems. 

Thesis Supervisor: Lawrence M. Lidsky 
Title: Professor of Nuclear Engineering 
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Chapter 1 


Introduction 


Inertial-electrostatic confinement (IEC) involves the creation of deep electrostatic poten- 
tial wells within a plasma in order to accelerate ions up to energies sufficient for fusion 
reactions to occur. These potential wells can be created and maintained by a slight ex- 
cess of electrons in a certain region of the plasma or by electrostatic grids. Typically 
such systems are arranged in a spherical geometry, as illustrated in Figure 1-1, taken 
from [1], Some of the earliest work on such systems was performed by Elmore, Tuck, and 
Watson [2], Farnsworth [3, 4], and Hirsch [5, 6]. A slightly different implemention has 
been described by Barnes, Nebel, and Turner[7]. 

Recently the fundamental IEC concept has been modified by Bussard [1, 8] to include 
a surrounding polyhedral cusp magnetic field in order to improve electron confinement; 
this type of system has been analyzed by Krall and Rosenberg [9, 10] and is depicted in 
Figure 1-2 (from ref. [1]). 

Another recent suggestion by Bussard [11], as well as Barnes and Turner [12], is 
to use driven acoustic standing waves to increase the average density in the core of the 
device. Such a process is illustrated in Figure 1-3 as originally presented in [11]. Following 
the convention of Bussard, this technique will be referred to as the inertial- collisional 
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-200 keV (TYPICAL) 


Inertial-electrostatic confinement: deep negative 
electric potential well 0 traps positive ion fuels (2) 
in spherical radial oscillations (3) until they make fu- 
sion reactions (4). 


Figure 1-1: Inertial-Electrostatic Confinement 
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ION 

SOURCE 



(D 


Inertial-electrostatic confinement: trapping well 
formed by energetic electron injection (T) into cusps 
of polyhedral magnetic fields (2) and ions fall into 
well and remain until reacted (3). 


Figure 1-2: Addition of Polyhedral Cusp Magnetic Field 
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compression effect, or the ICC effect. 


It has been suggested [1, 10] that IEC can maintain non-Maxwellian ion distributions 
at fusion reactor parameters; it has been further suggested in [1] that the two ion species 
may be kept at significantly different energies. Both of these properties would confer 
the ability for the fusion device to exploit resonance peaks in fusion cross sections more 
fully than other systems can. Such a machine would then be highly suitable for use with 
advanced fuels like D- 3 He and p- n B. Reactions using these fuels have the advantage of 
producing virtually all of their power in the form of charged particles which can be directly 
converted to electricity at very high efficiencies, as depicted in Figure 1-4 (from [1]). Since 
few neutrons are produced by burning these advanced fuels, radiation shielding require- 
ments and the activation of structural materials become much less worrisome problems. 

In addition to being able to use advanced fuels efficiently, IEC-based reactors might 
be able to offer higher power densities and simpler engineering designs than other fu- 
sion approaches. If these many advantages could be realized, IEC would overcome the 
objections that have been raised concerning more conventional fusion schemes [13]. 

The object of this paper is to examine various critical physics issues in as general 
a fashion as possible, so that the results will apply to a wide range of IEC systems 
and related variants. In particular, the potential problems which are analyzed include 
ion thermalization, ion losses, electron losses, bremsstrahlung emission, and synchrotron 
radiation losses. 

The following conventions are adopted throughout the paper, except where it is ex- 
plicitly stated otherwise. Temperatures and energies are in eV and all other quantities 
are in cgs units. If a species j is monoenergetic then Tj = \Ej. 
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-200 keV (TYPICAL) 

Direct electric conversion: reaction products CD are energetic charged particles, which escape against spherically sym- 
metric radial voltage gradient C|) to yield radiation-free direct electric power output (J). 


Figure 1-4: Direct Electric Conversion 
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Chapter 2 


General Physics Issues 


All of the issues examined in this section are effects that are essentially independent of 
the densities and density profiles which are used in the inertial-electrostatic fusion device 
(except for the very weak dependence on density contained in the Coulomb logarithm). 

In performing these generalized calculations, the following assumptions have been 
made: 

• There is no spatial variation of temperature or energy in the region of the device 
where the density is appreciable. (The potential well is assumed to have a broad, 
flat central region [1, 9].) 

• Other regions have low enough densities that most collision-related effects (eg. fu- 
sion, bremsstrahlung, and ion-electron heating) are negligible there. 

• Velocity distributions in the region of significant density fire approximately isotropic. 

• The spherical geometry returns scattered ions to the dense center of the device, 
so that only the parallel component of velocity-space diffusion must be considered. 
(In reality there will be some buildup of angular momentum due to perpendicular 
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velocity space diffusion in regions other than the exact center; this problem is treated 
in [10].) 

• There is no spatial variation in fuel stoichiometry (n,i/n; 2 , the ratio of the densities 
of the two ion species il and i2) in the region of appreciable density. 

• Quasineutrality (n e = Z^na + Z^rii 2 , in which n e is the electron density, while Zn 
and Z {2 are the charges of the two ion species) holds in the region of significant 
density. 


2.1 Fusion Power Density 


The gross power per volume produced by the fusion of two different ion species il and i2 
is 

= 1.602 - 10- 19 < (tv > Qm ini2 , (2.1) 

V cm 3 

where <7 is the fusion cross section, v is the characteristic velocity, Q is the energy (in 
eV) released per reaction, and n,i and n{ 2 are the densities of the two ion species. All 
quantities other than energy are in cgs units. In the case of fusion of a single ion species 
(eg. D-D reactions), tin n{ 2 in the above formula is replaced by | nf to avoid counting the 
same reactions twice. 


The fusion power can be expressed in terms of the ratio of the ion densities, x = 
«»i/n» 2 - If the first ion species has a charge of one and the second ion species has charge 
Z 2 , quasineutrality stipulates that n e = nji + Z 2 n{ 2 . Using these facts the fusion power 
density becomes 


P fus 

T~ 


1.602- 10 -19 < <tv > Q 


(z + z 2 y 


nz 


Watts 

cm 3 


( 2 . 2 ) 


Note that the fusion power is maximized for x = Z 2 . Contrary to what one may 
initially think, the power is not maximum for x = 1 since it is the total charge, not total 
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number of ions, which is being held constant as the fuel mixture is changed. The total 
charge is limited in general by the structure and strength of the confining electric and 
magnetic fields. 

For the purpose of comparisons with other characteristic times in the device, the 
characteristic fusion time of a test il ion with a member of the i2 ion species is readily 
defined (here Ui e // is the effective ion density seen by the test ion as it transits the 
system): 


Tfus = 


«»2 eff < (TV > 


(2.3) 


2.2 Energy Equilibration Between Ion Species 

It is worthwhile to check whether one ion species can be maintained at a significantly 
lower energy than the other ion species. To do so, it will be assumed that the il species 
is more energetic than the i2 species and that (at least to a first approximation) the 
standard Spitzer-type expression for interspecies energy transfer may be applied to this 
problem. 

As discussed in Glasstone and Lovberg [14] and Book [15], the heating of the i2 species 
by the il ions may be described by 


8 fcZ 4 %/ Zh Z?2 n t'l ^ A«'l— i2 frr m ^ 

~df = T 2x e (mara + maavTjV*" (T “ “ Tiz) ’ 


(2.4) 


in which all quantities are in cgs units. 


Converting all temperatures and energies to eV and evaluating the constants, one 
obtains: 


dT i2 _ , „ E , n-19 Vrnum^Z^ Zf 2 nn In A tl _ a 


dt 


= 1.75-10' 


{muTa + m i2 Ti x ) 3 / 2 


(T*-T i2 ). 


(2.5) 
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Because of the presumed Maxwellian distribution of i2 ions, the power density trans- 
ferred to them by the il ions is 


-P»l-«2 

V 


3 1 1 j-| — 19 \/m j i TTi^ 2 Zfi Zf 2 nn in A,i_i2 , 


( 2 . 6 ) 


It is now possible to consider two distinct cases. In the first case, T& is determined by 
balancing this heat transfer rate from il ions with the cooling effect due to the replacement 
of fused i2 ions with cold i2 ions. The second case is the situation in which the i2 ions axe 
somehow actively refrigerated to insure that they remain at very low energies. 

Proceeding with the evaluation of the first case, the cooling rate of i2 ions due to the 
replacement of fused ions is: 


P, tool 3^, 

~ y ~ ~ 2 T * 2 n * in *2 < > 


(2.7) 


The equilibrium temperature of the i2 species is determined by setting the total 
amounts of heating and cooling equal to each other. Since both the heating and cooling 
expressions have the same dependence on the ion densities, integrating them over the 
region of interest has no effect on the ratio between them. Solving for the equilibrium 
temperature and defining the ion mass as a multiple of the proton mass m p , mj = mm p , 
produces the result: 


T i2 = rn 


\ { 7.40 ■ 10 6 < av > (fiiiT i2 + (igTu) 3 / 2 
\J V-i\V-iiZ\ Z\ in Aa-,2 


( 2 . 8 ) 


For the fuels of interest (D-T, D- 3 He, p- xl B, etc.) typical values of the parameters 
are < av >~ 10 -16 - 10 -15 cm 3 /sec, Tn ~ 5 • 10 4 - 6 • 10 5 eV, and In A^_i2 ~ 15 — 20. 
By inserting any values within this range of parameters into the above formula, it is clear 
that 
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.957i! < Ti 2 < Tn • 


(2.9) 


Therefore from the evaluation of this first case it does not appear possible to keep 
one ion species at a significantly lower temperature or energy than the other without 
providing additional means of cooling the i2 species. 

Moving on to the second case in which a large temperature difference between the 
ion species is maintained by somehow actively refrigerating one species, it can be shown 
that the energy transfer rate required to sustain the nonequilibrium state would be pro- 
hibitively large. For this purpose assume that Tn > Ti 2 , so that the collisions between 
the two ion species occur at a velocity v ~ ^/3Tji/m,i. Coulomb collisions will then 
transfer energy between the species at the rate calculated above. Dividing this energy 
transfer rate by the fusion power, one obtains: 


i 20 • 10 -13 ( i2 
Pfus ' \m i2 / crQTa 


( 2 . 10 ) 


For a numerical estimate it is illustrative to use the case of p- n B reactions, for which 
it would be desirable to have high-energy protons (il species) and low-energy boron ions 
(i2 species). The peak of the fusion cross section, a a* 8 • 10 -2S cm 2 , occurs for a proton 
energy of about 620 keV, or Tn = |(620,000 eV). Estimating the Coulomb logarithm as 
approximately 15, the power ratio is found to be: 


P «l-» 2 ~ 
Pfus ~ 


1.4. 


( 2 . 11 ) 


If it were possible to operate the reactor with boron ions maintained at very low 
energies, the boron ions would siphon off more power froi the energetic protons than 
would be produced from the nuclear reactions. In order to keep the system operating 
precisely at the resonance peak of the reaction cross section, it would then be necessary 
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to renew the protons’ energy and refrigerate the boron ions; otherwise, the temperatures 
of the two species would rapidly equilibrate as in the first case. The task of continually 
adding so much energy to one species, subtracting it from the other, and extracting the 
resulting entropy from the system appears daunting at best. 

Thus both from the derivation of the natural energy equilibrium between ion species 
and from the calculated energy transfer required to keep the system in the nonequilibrium 
energy state, it does not appear to be possible to maintain one ion species at a significantly 
higher energy than the other. 

Even if both species can still be kept monoenergetic (but at approximately the same 
energy), the fact that < av > must be averaged over all collision angles then implies that 
it is impossible to exploit the resonance peaks of fusion cross sections (eg. the sharp peaks 
in the p- n B cross section) as fully as might be hoped. 


2.3 Ion Thermalization 


The ions will begin to evolve from their assumed initial monoenergetic distribution toward 
a Maxwellian distribution on a timescale characterized by the ion-ion collision time [16]: 


3A/B5T 1*1 , , VEElnl 

11 11 8x2- l 1 e 4 n i i e //lnA il -ii ' Z&nn *// In Aji-ii ’ 
The collision time may be compared with the fusion time, such that 


( 2 . 12 ) 


r «i-«i _ j 4 . 10 7 V^T Tf/ 2 <( tv > n i2 ' 

Tfu, ' Zfx InAii-ii n,i ’ 


(2.13) 


Once again using the typical values of < av 10~ 16 — 10~ 15 cm 3 /sec, Ti i rv 5 * 10 4 - 
6 • 10 5 eV, and In A tl _;i ~ 15 - 20, one finds that 
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~ 10 -3 - 10 -2 . 


(2.14) 


T fus 

Thus in the vicinity of the initial ion velocity, the ion distribution will begin to assume 
a Maxwellian form on a timescale which is two to three orders of magnitude faster than 
the fusion time. Of course the high-energy tail will require severed collision times to fill 
in, but even so it is apparent that the ion distributions will be essentially Maxwellian. It 
is also worth noting that the distribution will be truncated at the well depth energy, since 
the energetic tail can escape from the confining potential well. 

Because the ion distributions are Maxwellian to a good approximation, all of the 
fusion cross sections must be Maxwellian-averaged. One must therefore make do with the 
same < crv > values as are used in other fusion devices, and resonance peaks in the cross 
sections cannot be utilized more efficiently than in other types of reactors. In fact, because 
the high-energy tail of the Maxwellian is truncated at the well depth, the average fusion 
reactivities will actually be somewhat lower than truly Maxwellian-averaged quantities. 


2.4 Ion Upscattering Losses 

The ion losses due to radial energy upscattering can now be calculated. Ions can be lost 
either by completely escaping from the system or by climbing high enough in the well 
that the strong magnetic field near the plasma boundary deflects them into useless orbits. 
Both of these effects require that the ions be upscattered by a certain increment in energy. 

Consider the upscattering of a test il ion with charge Zn by field ions of both species. 
(The derivation that follows is an expanded version of an original derivation by M. Rosen- 
berg and N.A. Krall [17].) 

Sivukhin [18] gives suitable velocity-space diffusion coefficients for test ions (denoted 
by subscript t) amid a background of field ions which have an isotropic velocity distri- 
bution. The diffusion coefficients for monoenergetic background ions and Maxwellian 


19 



background ions are the same to within a factor of two; to emphasize the upscattering 
losses that can occur even if one assumes an initially monoenergetic distribution of ions, 
the coefficient appropriate to monoenergetic ions will be used here: 


47r2’, ? 1 e 4 lnA« ^ >72 2 

"II - 3m 2 iW 3 t Zj z i n i *ff V i 

■ < 2 -> 

in which it has been assumed that the two ion species have the same temperature. 

The diffusion coefficient may be related to Miyamoto’s collision time for il-il collisions: 



1 




(2.16) 


Using the velocity-space diffusion equation and recognizing that the test ion’s velocity 
can be related to the field ion velocity which it originally had by vn t = v,i(l + At>ii t/ v n)> 
one can calculate the time required for ions to be upscattered to velocities high enough 
for loss from the system: 


( Z i2 \ 2 n i2 mii | 

\Zn) nn T7Zj2j Tll-,1 ' 

(2.17) 

The loss time can be rewritten in terms of the energy upscattering, A E/E = 2Av/v, 
so that 





_ 3 /A E t ' 
~ 2 V E a 




-1 


Til—H 


(2.18) 
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This loss time may be checked by comparing it with the results of an earlier estimate 
of how quickly the tail of the Maxwellian is filled in from an initially monoene rgetic 
distribution, as performed by MacDonald, Rosenbluth, and Chuck [19]: 


m ii v fi t 


T tail 


12irZ t ^e 4 nii e // In A; 


ii 3 ^ T 2 E c ) 




(2.19) 


In the case that A E/E is of order unity (the typical value expected in IEC devices) 
and there is no i2 species, these two different expressions for the loss time are seen to be 
comparable. 

The fraction of ions that fuse is just the ratio of the loss and fusion times: 

Jloss _ o i in* 7 ( ( 1 4. — [l i ( \ ^ n >2 ^t'l \/ Mil 1 'ii £ ^ ^»2 

rjus ~ ‘ ' V E 0 ) V 2 E 0 ) l + UiJ Z 4 InA* n ia * 

(2.20) 

In calculating the amount of energy upscattering which is necessary for loss from the 
well, one must account for the fact that the well depth energy appears deeper to ions with 
Zn > 1. The condition for loss is that E 0 + A E t = ZnE we a , or 


A Et _ „ E we u , 

~e7 ~ Zii_ ir ~ 1 

A more complete calculation should also include cooling of the fast ions due to elec- 
tron drag. As a rough estimate of the importance of electron drag relative to ion-ion 
upscattering, one may consider the ratio of the ion-electron collision time to the ion-ion 
collision time, as defined in [16]. Considering only a single ion species for simplicity, it 
may be seen that 



Hi 

tu 


/v 


T e \ 3/2 


21VftZi (£) 


( 2 . 22 ) 


For the fuels and temperatures characteristic of the proposed IEC systems, the ion- 
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electron collision time is typically at least an order of magnitude larger than the ion-ion 
collision time. Therefore it appears doubtful that electron drag effects will substantially 
reduce the ion upscattering losses. 

Another possible mechanism for reducing the upscattering losses is cooling of the ions 
by charge exchange with neutrals. Unfortunately, while charge exchange may prevent ions 
from escaping the system, the neutrals themselves are free to escape, potentially carrying a 
sizeable amount of energy out of the confinement system. Considerable technical problems 
may also result from the pumping requirements necessary to collect the large quantities 
of escaping neutrals or ions and replace them with fresh ions. Of course, the presence of 
appreciable charge exchange effects would ultimately serve to degrade the monoenergetic 
ion distributions even more quickly than has already been calculated. Because of all of 
these reasons, the introduction of neutrals into the problem is not a useful solution. 


2.5 Bremsstrahlung Radiation Losses 

Maxon [20] gives the nonrelativistic and extreme relativistic limits for electron-ion and 
electron-electron bremsstrahlung and interpolates between these two limits to obtain ap- 
proximate radiation rates in the intermediate regime. His results may be used as a 
guideline for the necessary corrections for mildly relativistic (T e ~ 100 keV) electrons 
(consult McNally [21] for a commentary on the empirical expression for bremsstrahlung 
in this regime). Adding together the expressions for ion-electron and electron-electron 
bremsstrahlung and defining 


^e// — 


Ei x + z\ 


n e x + Z 2 

produces the result that the toted bremsstrahlung power per volume is: 


(2.23) 
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ftr 

V 


- l 1 + • 7936 ^ ? + 1,874 (^) ] + • or 


(2.24) 


T = 1.69- , [l + • + 1.874 (£)] + Wa '* ! 


•y/2 m e c 2 J cm 3 

(2.25) 


Once again using the expression given by Glasstone and Lovberg [14] and Book [15] for 
energy transfer between plasma species, the heating of electrons by ions may be described 

by 


dT, 

dt 


= 1.75 -10" 19 2 


yJm e mjZ\n{ In A c j _ _ v 
(miT e + m e Tj) 3 / 2 1 i e) 


(2.26) 


Because of the Maxwellian distribution of electrons, the power density transferred to 
them by the ions is § n e times the heating rate of Equation 2.26. After converting the 
power to Watts and assuming that T e » ^T,-, the result is 


^ = 7.61 • 10 -28 n* Y, - T.) , 

V i HiTe 2 cm3 


(2.27) 


in which the ion masses have been expressed in multiples of the proton mass, m; = Ai»7n p , 
temperature is still in eV, and density is in particles/cm 3 . 


As derived by Eosenbluth [23], there is a correction factor to this expression for ion- 
electron heating; the correction is caused by partial depletion of the electrons with veloc- 
ities smaller than the ion thermal velocity, with the net result that 


Pu 

V actual 


/27r 2 ^n,m e rA 2/3 ' 

V \3 5 / 4 n e rm T e ) 


( 2 . 28 ) 
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Furthermore, Dawson [22] notes that for relativistic electrons the ion-electron heating 
must be modified by a factor of (1 + 0.3 T e /m e c 2 ). After incorporating the corrections of 
Rosenbluth and Dawson, the heat transfer rate becomes 


Pie 

V 


7.61-10~ 28 n e 


V m e c 2 ) Y /tjTe /2 



(2.29) 


The equilibrium electron temperature is found by equating the power transferred 
to the electrons by ion-electron heating with the power lost by the electrons due to 
bremsstrahlung, synchrotron radiation, ion-electron cooling in the edge of the device, loss 
of electrons from the system, and other effects. The maximum possible bremsstrahlung 
rate may be obtained by neglecting all loss mechanisms except bremsstrahlung, thus 
producing the highest possible equilibrium electron temperature. In this approximation 
Pie = Pbr • (In the next chapter both synchrotron radiation and edge ion-electron heat 
transfer will be shown to be negligibly small compared with the bremsstrahlung and heat- 
ing effects considered here, so this approximation should come close to the actual answer.) 
Since the ion-electron heating and the bremsstahlung cooling both have the same depen- 
dence on the densities (with the exception of the Coulomb logarithm, which slowly varies 
from about 15 to approximately 20 over the range of the system), integrating over the 
system volume has essentially no effect. As a result, the equilibrium electron temperature 
(in eV) can be determined from the general equation: 


4.494. * fi + SS) V i [i - 

\ m e C 2 / Y Pi n e \ ^36 /*» n e Pf 


2/3 


= T 2 | 


z *ff 


1 + .7936-^- + 1.874 
m c c 2 


C Ti-T e ) 
3 T, 


(&)1 + £;&}• < 2 - 30 > 


After finding the equilibrium electron temperature, it may be used to calculate the 
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fraction of the gross fusion power output which is radiated away by bremsstrahlung: 


hr 

-P/uj 


1.06 -10" 13 


< <TV > 


+ l-874(^) 2 ] + ^^> 

^ (*+Z 2 ) J 


(2.31) 


As stated earlier, these values for the electron temperature and bremsstrahlung have 
been calculated assuming that bremsstrahlung is the dominant mechanism for cooling 
of the electrons. This assumption is justified because most of the other possible cooling 
effects (such as electron cusp losses) would introduce even greater power losses while 
cooling the electrons and reducing the bremsstrahlung power loss. Although cooling of 
hot electrons by cold ions in the edge would be a beneficial effect which would not cause 
further power losses, the very low densities in the edge region stipulate that the total 
electron cooling there will be much less than the total electron heating in the center of 
the device. 
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Chapter 3 


Design-Dependent Physics Issues 


Included in this section are effects that depend on the specific density, density profiles, 
and confinement system (eg. magnetic cusp, grids, etc.) which are employed. The first 
section will outline the specific assumptions made about the spatial density and energy 
profiles of the devices, and subsequent sections will use these profiles to calculate the 
magnitude of the various effects. 


3.1 Spatial Profiles 

8.1.1 Devices with Convergence-Limited Core Densities 

In the simplest EEC concepts, the core density is determined solely by the convergence of 
the spherical flow in the potential well. In the following calculations, it will be assumed 
that the device employs a single potential well. The theoretical analysis may be simplified 
by dividing the interior of the machine into three regions: the core (0 < r < r c ), the mantle 
(r c < r < r e ), and the edge (r e < r < R). Typically R ss 100r c and r e ~ 50 - 80r c . The 
following approximate forms for the particle densities and energies axe assumed. 
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Both the electron and ion densities are constant in the core, then because of conser- 
vation of particles in the nearly flat part of the potential well, they drop off like 1/r 2 in 
the mantle, and they finally reach a constant value in the edge: 


n e>i = < 


(n c ) e ,i(r c /r ) 2 

(n c ) e ,i(r c /r e ) 2 


0 < r < r c 
r c < r < r e 
r e <r < R. 


A graph of the density is given in Figure 3-1, taken from [9j. 


(3.1) 


Electrons in the core and mantle are heated enough by Coulomb friction with the 
energetic ions that they will tend to form a Maxwellian distribution with a temperature of 
T eo . As the electrons travel from the center to the edge of the well, they acquire additional 
energy corresponding to the well depth, so that the electron energy distribution is given 

by 


E e 


Zrp 
2 1 W 


0 < r < r e 


\T eo + E w f(r ) r e <r < R , 


(3.2) 


for which /(r) is some rapidly increasing function of r such that f(r e ) = 0 and f(R) = 1. 


Similarly, the ion energies have the following spatial variation: 



Zrp 

2*10 

§2»o - E w Zif(r ) 


0 < t < T e 

r e < r . 


(3.3) 


Figure 3-1 (from [9]) presents a graph of the potential well shape. 

To a first approximation bremsstrahlung, fusion, and ion-electron heating in the edge 
may be neglected because of the low densities and ion energies there. Using the fact that 
r c > r c , the following useful integral over the core and mantle regions is found: 



n 2 4irr 2 dr « — mi 2 r 2 . 
3 


(3.4) 
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3.1.2 Devices with Enhanced Core Densities 


Another class of EEC schemes is centered around increasing the core density beyond 
the convergence-limited value. References [11] and [12] propose to achieve this increase 
using acoustic standing waves; this phenomenon has been called the inertial- coliisional 
compression (ICC) effect. If the only role of the standing waves is to increase the core 
density or change the density profile, there will be no significant impact on the problems 
which have already been analyzed. The primary usefulness of the ICC effect lies being 
able to affect the total fusion power and hence perhaps the relative severity of electron 
cusp losses. 

If the root mean squared core density of an ICC device is enhanced beyond that of a 
non-ICC device by a factor of then the density profile may be approximated as 


n. 


(n c )e,i 0 < r < r c 

< ( n c)e,i{ r c/ r ) 2 /Z r e <r <r e 
. (n c ) eti (r c /r e ) 2 /Z r e < r < R . 


(3.5) 


For simplicity the energy profiles will be considered to remain approximately like those 
in purely convergence-limited machines. 

If the core density is significantly enhanced (£ > 1) via the ICC effect or other 
mechanisms, then essentially only the core will contribute to processes like fusion and 
bremsstrahlung; the rate of these processes in the mantle will be negligible by comparison. 
The net result of this fact is that 



n 2 4icr 2 dr 



2_3 

T 

c 1 c 


(3.6) 
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3.2 Relative Importance of Edge and Central Plasma Re' 
gions 


One of the key assumptions on which this entire analysis is based is that only the dense 
central region of the plasma contributes significantly to processes such as bremsstrahlung 
and ion-electron heat transfer; the edge region is assumed to make a negligible contribution 
to these processes. This assumption can now be justified. 

Effects such as bremsstrahlung and ion-electron heat transfer can be expressed as 
toteil powers and are proportional to the density squared times the volume of the region 
concerned. For a given process, the ratio of the powers of that process in the edge and 
central (combined core and mantle) regions is 

P'd ge ~ f*n 2 4Trr 2 dr 

Pcore+mantle So * ^4lCr 2 dr 

For a device with convergence-limited core density, Equations 3.1 and 3.4 may be used 
to express the ratio of integrals from Equation 3.7 as 


/JnW* 1 /r c \ ' 

fi‘nHirr*dr 4 \rj [\rj 


(3.8) 


For typical values of r e = 50r c and R = 2 r e , the ratio of edge effects to central region 
effects is 


n 2 4?rr 2 dr 

= °- 035 


(3.9) 


Of course, in computing the ratio of edge effects to central region effects for particular 
quantities like bremsstrahlung power or power transferred between ions and electrons, 
the exact answer will involve other numerical factors to account for parameters such as 
particle temperatures in the edge versus in the core, Coulomb logarithms in the two 
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different regions, etc. Precise evaluation of these factors requires detailed spatial profiles 
of the electron temperature and ion temperature in the edge region, but in general the 
net result of the factors will be to change the result above by at most about a factor of 
two or three. In conclusion, one finds that 


Pfdge 


Pi core+mantle 


~ 1(T 2 - 10 


-l 


(3.10) 


If the ICC effect is used to enhance the core density relative to the edge density, the 
ratio of the powers will be even smellier. It is clearly evident that for the purposes of the 
calculations presented in this paper, the edge region can be neglected in comparison with 
the centred region of the plasma for both non-ICC and ICC designs. 


3.3 Total Fusion Power 

Using the integred in Equation 3.4 for a machine with convergence-limited core densities, 
the toted fusion power is found to be 

P fua = 2.68 • 10~ 18 <crv> Qj^~J^ n le r c Watts, (3.11) 

where x = n^/n^, Zh = 1, and Zn = Z 2 . The fusion power is meodmized for x = Z%. 
For a device employing the ICC effect, Equation 3.11 is reduced by a factor of four: 

Pfu, = 6.70 • 10~ 19 <<7V> Q ^ n le r c Watts. (3.12) 
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3.4 Electron Cusp Losses 


An especially serious power loss mechanism is the loss of energetic electrons through the 
cusps of the confining magnetic field. At conditions of interest for a fusion reactor, the 
effective radius rjj of each cusp “hole” through which electrons escape is of the order of 
the electron gyroradius r e , so that tjj = knr e , where typically 1 < kg < 5 [24, 25]. The 
electron gyroradius, in turn, is given by the formula [15] 

r e = 2.38-^-p^- cm , (3.13) 

in which E e is in eV and B is in Gauss. Note that a factor of y/2 has been introduced into 
the usual formula because at the outer surface of the plasma the electrons are in directed 
motion [1]. 

Assuming a spherically symmetric distribution of electrons with essentially radial ve- 
locities, the number of electrons that escape each second will be \n e d ge e v e times the total 
area of the cusp holes. The factor of 1/2 is included because only half of the electrons are 
traveling outward. If each escaping electron carries away an amount of energy Ei oa , and 
there are N cusps (thus Nnrg is the total hole area), then the power loss due to escaping 
electrons will be 


P tlosa cusp — e'^eP'loaa • 


(3.14) 


Noting that Ei ots = E e = E we u, the well depth energy, and that v e = ^2E e /m e , and 
then expressing the power in Watts, energy in eV, and everything else in cgs units, one 
obtains 


P ! lost 


cusp 


1.69 . lo-n Mfagc&X* Watts , 
B £ 


(3.15) 
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If there is no ICC enhancement of the core density, then, the fraction of power which 
is lost because of escaping electrons is 


Peloas . ^ L = 6.30 * 


r i5/2 
"edge e& we U 


<orv> Q B 2 n\ e v\ 


(3.16) 


If the ICC effect is used to enhance the core density significantly, the vast majority 
of the fusion power will only come from the core, so the above expression should be 
multiplied by four. 

Now one malces the assumption that the edge of the plasma has /3 = 1 [26, 27], or 
B 2 = 8tt n e< ig e e Ewelh f° r Ewell hi ergs. Putting the energy in eV and substituting into 
Equation 3.16, the ratio of cusp loss power to fusion power can be simplified somewhat: 


Pe loss 


22 = 1 . 56 . 10 “#* 


i2 (* + Z2) 2 


H~ 


F 3/2 
A oell 


< <rv > Q nisi 


(3.17) 


One might think that using the ICC effect to increase the core density relative to the 
edge density would improve the ratio of cusp losses to fusion power, since the cusp losses 
occur at the edge and fusion occurs in or near the core. Yet as indicated in Equation 3.17, 
the constraint that /3 = 1 at the outer plasma boundary removed the dependence of the 
power loss fraction on the edge density, and thus on the density profile. It appears that 
using acoustic waves to alter the density profile of the device will create no significant 
improvement in reactor performance, provided that the waves only act to alter the density 
profile. The only critical parameter is the core density, which may be created via the ICC 
effect or simply by unaided ion flow convergence at the center of the device. 

The electron power loss fraction is minimized for x = at which 


P j? 3 / 2 

= 6.26- 10 18 JVJfe^ 2 , , . 

Pfu, H <av> Qn 2 e r 2 


(3.18) 


An expression for the characteristic electron loss time can also be derived. One begins 
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by noting that the fraction of electrons lost during each pass through the system will just 
be the total area of the cusp holes divided by the surface area of the machine. If this loss 
fraction is made small enough to be practical, then G e , the average number of transits 
an electron makes through the system before being lost, is well approximated as just the 
inverse of the loss fraction. Specifically, for N cusps: 


, 4 tR 2 R 2 B 2 

e ” Nvr 2 H ~ Nk 2 B E e • 


( 3 . 19 ) 


The loss time may be expressed as tj 044 = G e T tr , where r tr is the time required for a 
single transit through the system. This transit time is in turn r tr = 2 R/v av , where v av is 
the average electron velocity, or 


= A< 


m. 


( 3 . 20 ) 


with A some number of order unity to account for the faster electron speed near the edge 
of the well. 


Putting the above equations together, one obtains 


Te to., cusp ~ 1-68 • 10 " 8 N^AE^y/TTo ““ ’ (3 ' 21) 

where the temperature and energy are in eV and everything else is in cgs units. 


R 3 B 2 


One might be tempted to use electrostatic fields at the cusps in order to reduce the 
number of escaping electrons or the amount of energy which they carry away. Unfortu- 
nately, such techniques would also increase the ion losses or the energy carried away by 
each escaping ion, and so they are of little interest. 


A better method to reduce the electron power losses is to direct- convert the energy 
of escaping electrons into electricity. In particular, the addition of a “sideways” magnetic 
field outside each cusp would result in vxB forces which could separate outgoing escaping 
electrons from incoming fresh electrons. Then the outgoing electrons could be efficiently 
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directed around the electron guns, so that they will hit direct- converter grids and return 
most of their energy to the system. 


3.5 Electron Grid Losses 

Now consider the power losses that are caused by confining the particles with an elec- 
trostatic grid instead of a cusp. The ion losses on the grid can be minimized by making 
the grid bias large and positive, so that ions hitting the grid will possess essentially zero 
energy; one must then calculate only the electron grid losses. Assuming that the grid has 
radius r gT id and transparency to electrons rj ei and choosing the electron energy, velocity, 
and density to be evaluated at the grid, then one obtains: 


P c loss grid — (1 Ve^^^gridHgrid e^e gridEe grid • (3.22) 

Noting that v e = ^2E e fm e and expressing the power in Watts, energy in eV, and 
everything else in cgs units, one finds that 

Pe loss grid = 1.19 • l(T 10 (l - 7 u)r 2 grid n e g ridE Z J* rid Watts . (3.23) 

It is instructive to compare this expression for the grid losses with the earlier expression 
for total fusion power: 


3 /2 

P ! loss grid _ 4 45 . ^Q 7 C 31 jl 4 _ _ \ / T grid ^ 2 / n e grid \ grid 

Pfus * ' ‘ \ T c ) \ n ce ) < av > Qn ce r c ' 


(3.24) 


Taking r gr { d — R , ti c grid — He edge — 7i cc (r c /7' e )^, E c grid — Ewell i and x — Z 2 , this 
expression becomes 
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(3.25) 


Pe loss grid _ j , 
P fus 


10° Z- 


<‘-<y 




< (tv > Qn ce r e 


This ratio indicates that at typical reactor parameters the grid losses are several 
orders of magnitude greater than the fusion power. As an illustration of the optimum 
performance that can be expected, choosing Z 2 = 1, rj e = .99, R/r e = 2, E we u = 60,000 
eV, < av >= 10 -15 cm 3 /sec, Q = 2 • 10 7 eV, n ce = 10 18 cm -3 , and r c - 2 cm, one 
discovers that 


„ 3000 . 

Pf us 


(3.26) 


Not only are the electron losses tremendously greater than the fusion power, but one 
also has the inherent problem of cooling the grids. 

Although one might contemplate passing a current through the grid wires to create a 
magnetic field around them and reduce the number of particles striking the grid, this idea 
does not appear to be advisable. Magnetic fields strong enough to deflect particles from 
the grid wires would also interfere with the desired purely radial motion of the particles, 
thereby significantly reducing the degree of core convergence in the EEC system. 


3.6 Ion Grid Losses 

One could attempt to reduce the power losses by putting a large negative bias on the grids; 
then it will be the ions and not the electrons which constitute most of the power loss upon 
impact with the grids. By analogy with the electron calculation, it is straightforward 
to derive the ion losses caused by an electrostatic grid in the system. If the grid ion 
transparency is rj{ then 


Pi loss grid — (1 Vi)^^ r grid n i grid^i gridEi grid 


(3.27) 
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Noting that v,- = y/2EiJrni and expressing the power in Watts, energy in eV, and 
everything else in cgs units, one finds that 


Pi loss grid — 2.79 • 10 12 (1 — Tji)^ridP‘i gridE{ Igridl WattS . (3.28) 

For fii = 2 and all other parameters as before, the ion losses will be about 60 times 
smedler than the corresponding electron losses calculated above; the reason is simply that 
the ions are moving much more slowly than electrons of the same energy. Unfortunately, 
the grid losses are still much greater them the fusion power: 


Pj loss grid 
Pfus 

Because of the overwhelming power losses and cooling problems associated with grids, 
it appears preferable to use a different confinement technique such as the magnetic cusp 
system. 
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(3.29) 


3.7 Electron Thermalization 


Equation 3.21, the electron loss time for a cusp device with convergence-limited core 
density may be compared with the electron-electron collision time from [16]: 


. 25.8>i7Te 3/2 1q5 r e 3/2 

ee 16ir 3 / 2 e 4 n ee //lnA ee n ee ff In A ee ’ 


(3.30) 


where n ee // is the root mean square density experienced by an electron circulating 
through the plasma. If the density profile may be approximated by Equation 3.1 with 
R « 100r c , n e ff will be about an order of magnitude smaller than the core density. 


It can be determined whether the electrons will be significantly thermalized by con- 
sidering the ratio of the two times, 
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(3.31) 


= , 9 . 10 u Tjltl&AZ* l 

Te loss cusp I® Ajjflj e ffB?B 2 

For the typical parameters N = 8, kjj = 2, A ~ 1.5, T e< > ~ 2 • 10 4 - 1.5 • 10 s eV, 
E wM ~ 6 • 10 4 - 9 • 10 5 eV, In A ee ~ 15 - 20, n e e// ~ 10 16 - 10 17 cm" 3 , iZ = 150 - 300 
cm, and B ~ 10 4 - 10 6 Gauss, the ratio is found to be in the range 

— — ~ 10 -6 - 10 -3 . (3.32) 

T e loss cusp 

Thus it is readily apparent that electrons in the center of the EEC device will form an 
essentially Maxwellian distribution. 


3.8 Synchrotron Radiation Losses 


In calculating the electron temperature and bremsstrahlung losses in the previous chapter, 
the effects of synchrotron radiation were assumed to be negligibly small. This assumption 
will now be justified. 

The power density of emitted synchrotron radiation is given in [16] as: 


Psyn iMPnc ( T e \ 5 f T e Y 

V 3 m^c 3 \m e c 2 ) . 2 \m e c 2 /. 


(3.33) 


Evaluating the constants, defining V 3yn to be the plasma volume which is under the 
influence of the magnetic field and emitting synchrotron radiation, and letting / represent 
the fraction of the radiation which is actually lost (not reflected back into the plasma and 
reabsorbed there), the synchrotron power becomes 


P, yn = 6.21 • 10 -28 J3 2 n e T„ 


1 + 


2 \m e c 2 J. 


fV 3 ^ Watts. 


(3.34) 
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In a diamagnetic EEC plasma, synchrotron radiation will only come from the outer 
layer of the plasma. Electron diamagnetism prevents the external magnetic field from 
penetrating more than a few electron gyroradii into the plasma [26]. Using the fact 
that in this outer layer B 2 /8ir = n e T e and defining the layer’s thickness to be kur e , the 
synchrotron power is found to be 


JV = 2.50< 


10 ~ 38 fTZ edge 


1 + 


!( : 


— — y-')] n; 
m e c 2 )\ 


Idge e 4*R 2 k H r e Watts. 


(3.35) 


The condition that B 2 /8w = n e T e allows the electron gyroradius to be rewritten in 
terms of the density: 


r e ~ 2.38 


•yZ-Tg ~tdgt 


B 


cm = 3.75- 10 s - 


yj^edge * 


cm. 


(3.36) 


With the aid of the relations {n e d ge e /n ce ) = ( r c /r e ) 2 and T ee d ge = \E W eiu the ratio 
of the total synchrotron power to the total bremsstrahlung power may be estimated: 


ayn 


0.19/te (^) (0 


R \ 2 Elil 


(3.37) 


Pbr \ r eJ \Te / ]j leo T c yjTl ce 

(For this estimate the relativistic corrections to the synchrotron and bremsstrahlung losses 
were neglected, since they were of the same order of magnitude.) 

Even at the rather extreme parameters of E we u = 10 6 eV, r c = 2 cm, n ee = 5 • 10 17 
cm -3 , E we u/T eo ~ 7, R/r e = 2, and r c /r e = 1/50, the ratio is only: 


•* syn 

Pbr 


0.03 fk H . 


(3.38) 


Since / < 1 and kjj should at most be 4 or 5, it is clear that bremstrahlung, not 
synchrotron radiation, will be the dominant radiation loss mechanism. (Note that there 
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is the luxury of further reducing the synchrotron losses by reflecting and reabsorbing most 
of the radiation in the plasma, thus making / much smellier than 1.) 

This analysis of the synchrotron radiation losses has assumed that the plasma diamag- 
netically excludes the magnetic field except in a very thin sheath at the plasma surface. 
Even if such diamagnetic effects do not occur, the synchrotron losses will be limited by 
the fact that the vacuum magnetic field of the cusp system varies as ( r/R ) n , where n > 3 
[10]. Thus a strong magnetic field will still exist only near the plasma surface, although 
it might penetrate far enough so that the sychrotron losses become comparable to the 
bremsstrahlung losses. (In that event, however, electron cusp losses would become so 
severe that they would be a much more pressing concern than synchrotron losses.) 
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Chapter 4 


Performance for Various Fuels 


The following tables present parameters and results for IEC reactors using various fuels. 
The reactors are assumed to use a magnetic cusp confinement system but not to utilize 
the ICC effect to enhance the core densities beyond those of normal convergence-limited 
flow. In each case the parameters have been chosen so that the reactor performance is 
approximately optimized. 

The ions are assumed to have an initially monoenergetic distribution at an energy 
Ei 0 = E we u/2. (Of course to produce such an initial distribution one would have to 
surmount the difficulties of very accurately injecting ions that deeply into the well, but 
for the present calculations those technical problems are neglected.) The ions will begin 
to evolve toward a Maxwellian distribution with T; 0 = (2/3 )Ei 0 on a timescale of r^. As 
is indicated for each fuel species in the tables, the thermalization of a test ion typically 
occurs two to three orders of magnitude more rapidly than the fusion of that test ion. 
Thus the ion distributions in the center of the device will be essentially Maxwellian, except 
that the high-energy tail will be truncated at the well depth. 

Neglecting this truncation, the appropriate reactivity to use in these calculations is the 
beam- Maxwellian quantity. However, under these conditions where the beam energy is 3/2 
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of the Maxwellian temperature, the reactivity can be well approximated to within a few 
percent by the Maxwellian-averaged reactivity; hence the Maxwellian-averaged < av > 
values are used for simplicity. Cross section data is drawn from references [28], [29], and 
[30]. 


Parameter 

D-T 

D- 3 He 

D-D 

Swell 

60 keV 

210 keV 

300 keV 

T io 

20 keV 

70 keV 

100 keV 

T'o 

18 keV 

56 keV 

76 keV 

n ce 

5 • 10 17 cm -3 

5 • 10 17 cm -3 

1 • 10 18 cm" 3 

B 

2.2 T 

4.1 T 

7.0 T 

Fuel mixture 

1:1 

1:1 

- 

r c 

1.5 cm 

2.5 cm 

2.5 cm 

(III A ) average 

16 

16 

16 

< (TV > f us 

(10 -16 cm 3 /s) 

4.24 

1.02 

0.45 

Q 

17.6 MeV 

18.3 MeV 

3.7 MeV 

Ncuapa 

8 

8 

8 

kH 

2 

2 

2 

P fua 

4.2 GW(t) 

2.2 GW(t) 

3.5 GW(t) 

Pneutrona / Pfua 

0.80 

0.01 

0.36 

Pbrem/P fua 

0.008 

0.24 

0.52 

Pe loss cusp/ P. fus 

(from Eq. A, 5) 

0.11 

1.32 

1.42 

Pe loss cusp/ P, fus 

(from Eq. 3.17) 

0.47 

5.76 

6.39 

loaa/'Tfua 

D: 5 • 10- 3 
T: 4 • 10~ 3 

D: 3 • 10" 3 
3 He: 3 • 10" 2 

CO 

o 

1 

T ii/ r fua 

D: 1 • 10~ 3 
T: 2 • 10 -3 

D: 2 • 10" 3 
3 He: 2 • 10~ 4 

2 • 10- 3 

'Tee/'Te loaa cuap 

1 • 10 -5 

3 • 10~ 5 

1 • lO" 5 


Table 4-1: IEC Reactors Utilizing Deuteron-Based Fuels 



Parameter 

p- n B 

p- 6 Li 

Ewell 

900 keV 

1.5 MeV 

T io 

300 keV 

500 keV 

Teo 

138 keV 

206 keV 

Tice 

5 • 10 17 cm -3 

kbebxsb 

B 

8.5 T 

15.5 T 

Fuel mixture 

5:1 p: ll B 

3:1 p: 6 Li 

r c 

3 cm 

3 cm 

isia—i 

16 

16 

< (TV > f us 

(10 -16 cm 3 /s) 

2.39 

1.1 

Q 

8.7 MeV 

4.0 MeV 

N 

■L'cuaps 

8 

8 


2 

2 

P' fus 

1.9 GW(t) 

2.7 GW(t) 


< 10" 3 

- 

■E3SE991I 

1.73 

5.36 

Be lose cusp / P fua 

(fromEq. A.5) 

10.6 

15.1 


61.1 

92.5 

loss /^f us 

p: 2 • 10“ 2 

U B: 6 


Tii/Tfu, 


■KBfiBfl 


1 • 10" 4 



Table 4-2: IEC Reactors Utilizing Proton-Based Fuels 


The tables also compare the time required for ions to be lost via upscattering into the 
high-energy tail with the time required for them to fuse. Since £ to = E we u/2, Z — 1 ions 
will be lost when A Et/E 0 = 1. With the exception of ions such as n B and 6 Li which see 
a well much deeper than twice their initial energy, the ions will escape from the system 
far more rapidly than they will fuse. Typically between 30 and 300 ions will escape for 
every ion that fuses. As has already been noted, even if charge exchange can solve the 
problem of ion losses, it will create the difficulty of the loss of fast neutrals and compound 
the problem of rapid ion thermalization. 
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Because both of the ion species are nearly Maxwellian with equal temperatures, the 
averaged fusion cross sections will be considerably smaller them they would be if the ions 
could be kept monoenergetic at the resonance peak energy. In fact, since the high-energy 
ion tail escapes from the system, the average reactivities of fuels in an IEC device will 
be somewhat less them those in a fusion reactor which can confine the hot ion tail of the 
Maxwellian. For simplicity, Maxwellian-averaged reactivies were used in calculating the 
fusion power, but it must be remembered that the true reactivity will be smaller and so 
the power loss fractions will be somewhat larger than shown. 

While the bremsstrahlung power loss is quite small for D-T and more or less tolerable 
for both D-D and D- 3 He, it is found to be prohibitive for the other fuels, since the high 
ion energies in the center of the device lead to high electron temperatures there as well. 
Because of the truncated Maxwellian ion distributions, the bremsstrahlung/fusion ratios 
will be roughly equal to or perhaps even worse than those of other reactors burning these 
fuels. 

The strength of the cusp magnetic field is calculated assuming that /? = 1 at the outer 
plasma surface [26, 27], or B 2 = 8rni e d ge e E we u with E we u in ergs. As given in Equation 
3.1, n e dgee = n ce (r c /r e ) 2 . For the purposes of these calculations, it has been assumed 
that r e = 50r c . 

There are two different ways to calculate the electron power losses. Equation 3.17 
gives more pessimistic results than Equation A.5, but the latter equation is based on 
many assumptions about the potential well and density profiles, whereas the former is 
not. The tables give the ratio of electron loss power to fusion power as determined by 
each method; in using Equation A.5, it was assumed that A — 1.5 and R = 2 r e . 

Even by using the more optimistic answer and adding direct converters with 50-60% 
efficiency, the electron losses appear to be intolerably large for fuels other than D-T. It 
should also be noted that this calculation was based on the optimistic assumptions that 
kjj could be kept as small as 2, the effective cusp number would not be larger than 8, 
very high core densities could be achieved, and Ohmic power losses in the field coils could 
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be neglected. Under actual conditions, the losses will probably be even more severe than 
those calculated here. 

(It would be possible to reduce the electron cusp losses if the outer layer of the diamag- 
netic plasma could be maintained in equilibrium with /3 < 1, so that higher magnetic field 
strengths could be used. However, the behavior of the outer sheath of the diamagnetic 
plasma is poorly understood, and the plasma might simply adjust itself to keep /3 = 1, 
as assumed in [26]. In any event, the magnetic field strengths indicated in the tables are 
already quite large, so it would be rather difficult to increase them much more.) 
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Chapter 5 


Conclusions 


The suitability of various implementations of inertial-electrostatic confinement (IEC) sys- 
tems for use as D-T, D-D, D- 3 He, p- n B, and p- 6 Li reactors has been examined. It has 
been shown that while an IEC reactor would have the advantages of high power densi- 
ties and relatively simple engineering design when compared with other fusion schemes, 
it suffers from several flaws. These problems are ion thermalization and upscattering, 
bremsstrahlung, and electron cusp losses. 


5.1 Ion Thermalization and Upscattering 

The problem of ion thermalization and upscattering can be described in a straightforward 
manner. A test ion is injected into the well at the desired energy and begins to oscillate 
through the dense core, out toward the plasma edge, and back again. Collisions with 
other ions, all presumably starting at the same energy, will cause the test ion to diffuse in 
velocity space. The perpendicular velocity-space diffusion is less important than parallel 
diffusion, since most scattering occurs near the device center, and the well then returns 
the scattered particles to the center for another try. (There will be some core spreading 
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due to collisions away from the exact center, but the calculations in [10] show that the 
spreading is tolerably small.) For this reason, the present analysis has focused on parallel 
velocity-space diffusion, or energy up- and downscattering. 

Since the test ion only spends a fraction of its time in the core, it would be incorrect 
to compute the ion-ion collision time using the core density; rather, one must use some 
sort of average density seen by the ion as it transits the entire system. Because the 
average density is typically much, smaller than the core density, the thermalization and 
upscattering times will be significantly lengthened. 

Next the characteristic thermalization and upscattering times for the test ion must be 
compared with the characteristic fusion time for that ion, in order to determine whether 
the ion is likely to fuse first or be upscattered out of the well first. Since the test ion 
spends only a small fraction of its time in the core, its fusion time should be lengthened by 
the same factor and for the same reasons as the collision time. When one takes the ratio 
of the fusion and thermalization times, which axe both inversely proportional to density, 
the two factors cancel each other. One finds the ratio to be the same as in other fusion 
reactor designs, namely that an ion thermalizes about two to three orders of magnitude 
faster them it fuses. Therefore, it appears that the ion velocity distribution in the broad 
flat bottom of the well will look essentially like a Maxwellian truncated at the well depth 
rather than the desired monoenergetic distribution. 

If the initial ion energy is not much smaller than the potential well depth, say half 
of the well depth, then ions will not have to be scattered terribly far out into the tail to 
be lost. They can cover this comparatively small distance in velocity space in just a few 
ion-ion collision times. Thus not only will the ions thermalize far more rapidly than they 
will fuse, but they will also escape the well much more rapidly than they fuse (except for 
ions like n B which see a much deeper well). 

Furthermore, energy is transferred between the two ion species on a time scale roughly 
comparable to the thermalization time of each individual species, so that it would not 
appear possible to maintain the two ion species at significantly different energies or tem- 
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peratures in order to take better advantage of the resonance peaks in the reaction cross 
sections. It was shown that even if the ion species could be kept at vastly different ener- 
gies, the heat transfer between them would drain energy from the more energetic species 
at a rate comparable to the total fusion power. This heat transfer rate would necessitate 
the introduction of some method for returning this large amount of power from the lower- 
energy ion species back to the higher-energy species and somehow pumping the newly 
generated entropy out of the system; in practice it would be nearly impossible to fulfill 
this requirement. 


5.2 Bremsstrahiung 

The radiated bremsstrahiung power will of course depend on the rate of the energy transfer 
between ions and electrons which occurs primarily near the dense core. Just as the ion-ion 
collision time is much shorter than the fusion time, it can also be easily shown that the 
electron-electron collision time is many orders of magnitude shorter than the electron loss 
time in the magnetic cusp confinement system. Therefore both the ions and electrons 
will have roughly Maxwellian distributions in the device center; the temperature of the 
electrons relative to that of the ions must then be determined. 

For a given ion temperature, the equilibrium value of the electron temperature is 
obtained by equating the standard Spitzer-type ion-electron heat transfer rate with the 
bremsstrahiung cooling rate of the electrons; this derivation yields electron temperatures 
which are typically at least half of the ion temperature. As a result, at the very high 
ion energies necessary for p- n B and p- 6 Li reactions, the bremsstrahiung is prohibitively 
large, just as it is in more conventional reactor designs. 

One possible solution to the bremsstrahiung problem for advanced fuels is that the 
ion-electron heat transfer rate for highly non-Maxwellian species or species at widely 
differing energies might be significantly lower than the standard Spitzer heat transfer 
rate. However, this solution would require methods of maintaining nonthermal particle 
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distributions, a feat which IEC does not presently appear able to accomplish. 

5.3 Electron Cusp Losses 

In the analysis of the power loss due to electrons leaking through the point cusps, many 
optimistic assumptions were made. It was assumed in the calculations that strong plasma 
diamagnetism restricts the magnetic field to a thin outer sheath of the plasma and only 
allows electrons to leak through cusp holes which have radii of a few electron gyroradii. 
Furthermore, the loss hole radius was chosen to be only twice the electron gyroradius 
(according to work by Grad [24] and Grossman [25] it could be as large as five). Another 
assumption was that there were only eight point cusps, even though the effective number 
of cusps may be significantly larger due to practical limitations on the magnetic field 
geometry, or it may be necessary to increase the number of cusps to keep the system 
more nearly isotropic. Using these assumptions an expression for the electron power loss 
was derived, and it was found to be roughly proportional to E^ u , where E we u is the well 
depth. With the appropriate choice of parameters, and provided that the above-stated 
optimistic assumptions hold, the cusp power losses for D-T are tolerable, because a well 
depth of only a few tens of keV is required. However, the significantly greater well depths 
required for all other fuels cause their electron power losses to be prohibitively large. 

Future EEC research should more closely examine the diamagnetic “whiffle belli” and 
sheath effects in magnetic cusp confinement systems, as these phenomena are poorly 
understood at the present time. If the outer plasma sheath can be kept at /3 < 1, the 
electron confinement could be somewhat better than predicted here. On the other hand, 
if the loss hole radius cannot be made as small as a couple of electron gyroradii, the 
electron confinement will be much worse than was calculated. 

If electrostatic grids are used instead of magnetic cusps, the electron losses should be 
orders of magnitude worse, large numbers of ions would also be lost by collisions with the 
grids, and the severe problem of grid heating would also arise. While grids are convenient 
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for small-scale experiments, they do not appear to be desirable in actual IEC reactors. 


5.4 Acoustic- Wave Compression of the Core 


Although the use of acoustic standing waves to increase the core density and/or alter the 
density profile has been proposed in both [11] and [12], it appears that such a phenomenon 
could do little to improve the fundamental problems noted above (and it may even have 
a detrimental impact). For example, the ratios of ion thermalization and upscattering 
times to the fusion time are independent of both the core density and the spatial profile 
of the density in the reactor, and so they would remain unaffected by the so-called ICC 
effect. Likewise the ratio of bremsstrahlung power to fusion power would also remain the 
same. 

As noted previously, one might think that using the ICC effect to increase the core 
density relative to the edge density would improve the ratio of cusp losses to fusion power, 
since the cusp losses occur at the edge and fusion occurs in or near the core. Yet as it was 
shown, the constraint that (5 = 1 at the outer plasma boundary removed the dependence 
of the power loss fraction on the edge density, and thus on the density profile. 

The only critical parameter is the core density, which may be created via the ICC 
effect or simply by unaided ion flow convergence at the center of the device. Obviously 
the primary effect of altering the core density will be to change the fusion power density 
and total fusion power. 

Even if it were quite desirable to employ the ICC effect, it is far from certain that 
the acoustic waves will work as expected to compress the core. If the ICC effect does 
indeed occur, it is highly questionable whether it can achieve the necessary many-fold 
compression without simultaneously degrading the central ion convergence and defeating 
the purpose of its use. 
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5.5 Other Potential Problems 


There are several other issues which were not examined in this paper but which would 
need to be carefully considered in future IEC work. These areas include determining 
the limitations on maximum core density and more closely scrutinizing the rate of core 
spreading due to angular momentum buildup. One would also have to perform analyses of 
counterstreaming and Weibel instabilities, taking into account the fundamental nonlinear, 
nonlocal nature of the problem. Another question is whether part or all of the potential 
well will eventually fill in due to background neutrals and other effects. In addition to these 
physics issues, there are serious technological problems which must be explored, such as 
finding suitable techniques for accurately fueling deep inside the well and designing direct 
converters appropriate for the spherical geometry of IEC devices. 


In conclusion, it is hoped that discussion of these apparent problems will in the future 
lead to the discovery of methods which can circumvent them, allowing IEC devices to 
maintain energetic non-Maxwellian ion populations with relatively cold electrons and to 
offer a good power balance even for advanced fuels. 


51 



Appendix A 

Alternate Derivation of Cusp 
Losses 


An alternate way to derive the power loss due to electrons leaking through the mag- 
netic cusps is to make use of tj OS3 , the characteristic loss time of the electrons, as given in 
Equation 3.21. If the toted electron population in the machine is N e and the energy per 
lost electron is E we u (in ergs), then the power loss due to escaping electrons is 


> . _ N e E we u 

e loss cusp — 

Tloss 


(A.1) 


Using the density profile of Equation 3.1, the total electron population is found to be 


N e = -irn ce rfr e 




(A.2) 


For typiced length ratios within the machine, r c /r e <C 1, so it may be neglected compared 
with the other two terms within the brackets. 
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Putting the above equations together and expressing E we u in eV, one obtains 


P , .... = 3.98 ■ I<r”j4 [2 + (0 3 | Watts. (A.3) 

Setting /3 = 1 at the plasma edge so that B 2 = 8irn e< i g e e E we n, and using n e <i 9 ee = 
n ce (r c /r e ) 2 , the cusp power loss becomes 


Pe loss cusp — 0.988A |l + 2 


fE we ii\/T, 'Z Watts. 


The ratio of the electron loss power to the fusion power of an IEC device with 
convergence-limited core densities is 


= 3.69 • 10 17 A [l + 2 f-jj-Vl ( z + Z *)* ^k 2 H E w ± uVT^ . 
Pfus [ \RJ J * <<TV> Qn^rf 


(A.5) 


Several factors in this expression are determined by the shape and depth of the electro- 
static well. Using typical values of Ak 1.5, r e /R ns 0.5, and T eo ns E we ii/8, Equation A.5 
reduces to 


Pe lota cusp ^ ^ < ^Q 17 ^ 

Pf us ~ x <crv> Qn\ t r\ ' 

This answer is about 4 times smaller than the result of the electron loss derivation 
that was presented earlier. However, it is dependent on the exact potential well profile, 
density profile, and relative temperatures of the different particle species, whereas the 
earlier result was not. In any event, there is enough uncertaintly in kg, which enters 
both equations as a square, that a factor of 4 is comparatively unimportant. Whichever 
derivation one uses, the electron losses are sizable for D-T and intolerable for D- 3 He and 
p- 11 B, even when one optimistically assumes some relatively small value for kg. 
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